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ABSTRACT 

We introduce a new technique to detect the discrete temperature steps that cosmic 
strings might have left in the cosmic microwave background (CMB) anisotropy map. 
The technique provides a validity test on the pattern search of cosmic strings that 
could serve as the groundwork for future pattern searches. The detecting power of the 
technique is only constrained by two unavoidable features of CMB data: (1) the finite 
pixelization of the sky map and (2) the Gaussian fluctuation from instrumental noise 
and primordial anisotropy. We set the upper limit on the cosmic string parameter as 
G\i < 3.7 x 10~ 6 at the 95% confidence level (CL) and find that the amplitude of the 
temperature step has to be greater than 44/xK in order to be detected for the Wilkinson 
Microwave Anisotropy Probe (WMAP) 3 year data. 

Subject headings: CMB anisotropy, cosmic strings 



1. Introduction and Modeling 



Cosmic string is one of the relic structures that are predicted to be produced in the course 
of symmetry breaking in the hot, early universe, whose discovery will probably be an impor- 
tant landmark for the high energy physics. The quest for cosmic strings has been conducted 
in two ways: statistical method and direct search for individual cosmic string. On the statis- 
tical side, many studies agreed that the contributi on from cosmic strings to stat i stical observ- 
ables such as the power spectrum is at most 10% (jPogosian et al. 20031 : IWu 20051 : iFraisse 20051 : 
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square patch of constant background Gaussian fluctuation temperature step 
CMB sky map temperature (CMB + noise) by cosmic string 



Fig. 1. — Decomposition of a CMB anisotropy sky map patch. A square patch is decomposed into 
the uniform background temperature (To) + Gaussian fluctuation (o"g)+ discrete temperature step 
(±A). The number of pixels in a patch is N 2 . 



Davis k Kibble 20051 : IWyman et al. 20051 ). reconfirming that cosmic strings played a minor role, if 
any, in making the universe. Other workers have set upper limits on the cosmic string parameter 



Gu (IHindmarsh 1994 ; iPerivolaropoulos fc Simatos 1998 ; IJeong k Smoot 20051 ; iLo k Wright 20051 ; 
Seljak fc Slosar 20061 : iPogosian et al. 20061 : iFraisse 20061 ). In this Letter, we introduce a technique 
for estimating how strong a signal from cosmic strings has to be in order to be identified unequiv- 
ocally. The detecting power of this technique is calibrated by applying it to simulated anisotropy 
maps based on reasonable modeling. In the last part of this Letter, we apply this technique to 
the WMAP 3 year W-band data selo and analyze the implication of the results. A cosmic string 
can leave discrete tempe rature steps in a CMB anisotropy map due to the Kaiser-Stebbins effect 
(IKaiser k Stebbins 1984 ) , with the height of a step 5T given by 



5T = 8TTGn-f s p s Th -(vxs) 



(1) 



where T ~ 2.7251T is the universal background temperature of CMB. However, those step patterns 
are probably obscured by instrumental noises and other physical structures of anisotropies. We 
consider a square patch of the CMB anisotropy sky map that extends to a size of the horizon at 
the time of recombination so that the conic spacetime formed by a cosmic string can be in effect 
within the region. Then another square patch of the CMB anisotropy sky map that contains a 
segment of moving cosmic string can be decomposed into three parts, as illustrated in Figure [TJ 
(1) the uniform background temperature (To), (2) the Gaussian fluctuation (variance ctq), and 
(3) the discrete temperature step (±A) from a moving cosmic string. The uniform background 
temperature To comes from superhorizon-scale primordial fluctuations. Since we pick a horizon- 
sized region, the super-horizon fluctuations will appear to be a constant temperature shift for the 
whole patch. Smaller scale (sub-horizon) fluctuations and instrumental noise add up incoherently 



to form a Gaussian fluctuation with the variance Oq 



'CMB 



<oisei where °CMB and ° noise are 



1 See http://lambda.gsfc.nasa.gov/product/map/ 
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variances for fluctuations of primordial origin and instrumental noise, respectively. We introduce 
five parameters that characterize a square patch of a string-embedded sky map, To, era, A, p 
(blueshifted pixels/total pixels), and 9 (orientation of step). We concoct a simulated patch of the 
CMB anisotropy map by assigning arbitrary legitimate values to the parameters and adding the 
three components illustrated in Figure [TJ To recover these parameters from the CMB anisotropy 
map where the step pattern is intermixed, we employ five observables (p, A, 6, uq, To) that can 
be expressed in terms of 



mean 



T = N- 2 Y,Ti (2) 

% 

dipole moment : d = iV~ 3 T%£% (3) 

i 

A -inertia : A = 4iV~ 3 ^T |d • r*| / |d| (4) 

i 

variance : a 2 = N~ 2 ^ (T - r) 2 (5) 

i 

where T represents the signal at the iih pixel and r$ is the two-dimensional gridded position vector 
for the ith pixel defined in a square patch. In the absence of the background fluctuation (<tg = 0), 
the following analytic expressions for the characteristic parameters return the exact input values: 

2d- i(|A-r| + A-r) 
P = 4d-|A-r| (6) 

A = (' A - T '- 4d ) 2 (7) 
2(2d-|A-r|) K> 

9 = tan^ 1 (d y /d x ) (8) 
a G = {a 2 -4A 2 p(l-p)) 1/2 (9) 
T = r-A(2p-l). (10) 

With the background fluctuation turned on, an observed characteristic parameter Qj is expressed 
as an unbiased estimate Qi plus a Gaussian error ag_ i (Qi = p, A, 9, aa, To), 

Qi,obs = Qi±<r Qi . (11) 

The errors (<r p , a a, &6i a a a , ct ) present in the equations (fTTI) are measures of how precisely the 
information on a temperature step screened by background fluctuation is recovered. Figure [5] shows 
the behaviors of errors as functions of A for three different pixel numbers for a patch. The standard 
deviation of the background fluctuation chosen here is the mean standard deviation of patches with 
angular radius Or = 1.8° (the angular size of the horizon at the recombination for the ACDM 
model) for the WMAP 3 year W-band data set with the Kp2 mask applied. The most pronounced 
feature of the graphs is that the error bars are wildly undulating for small A compared to ac, with 
attenuating envelopes as A increases. The algorithm works poorly for low A compared to ctg, such 
as the ±90° error on the orientation estimate. Even the breakdowns of the algorithm do happen 
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Parameter Input Output (1 a) 
Patch with 8x8 pixels 

p 05 0.5 ± 0.1 

A (pK) 80.0 83.4 ± 16.6 

B (deg) 0.0 -0.1 ± 12.3 

CT G (/iK) 116.0 113.0 ± 12.0 

T (/xK) 10.0 9.7 ± 18.9 

Patch with 14 x 14 pixels 

p 05 0.5 ± 0.1 

A (/iK) 45.0 47.1 ± 9.5 

6 (deg) 0.0 -0.2 ± 12.6 

a G (fjK) 116.0 115.0 ± 6.2 

Tp(fiK) 1O0 9.6 ± 10.9 

Patch with 28 x 28 pixels 

p 05 0.5 ± 0.1 

A (/iK) 25.0 25.8 ± 4.7 

(deg) 0.0 0.2 ± 11.2 

CT G (/iK) 116.0 115.7 ± 3.0 

T ([iK) 10.0 9.9 ± 5.3 



Table 1: Simulation results at the critical values of A above which cjas are small enough and do 
not decrease any more. 

in the shaded region, resulting in erroneous parameter values such as negative A or p not in the 
range between and 1. Simulations with collapsed results are not included in the statistics shown in 
Figure[2]since those cases are evidently the ones with insufficiently strong signals against background 
fluctuation. As A increases, we begin obtaining the computed parameters that are very close to 
the true values with acceptable errors (dubbed the "good" results), and it allows us to recover 
the temperature step parameters faithfully. We also find from Figure [2] that a patch with more 
pixels works better with a faster attenuation of uncertainties. Table [T] displays the performances 
of the algorithm for the values of A at the borders, above which errors for the temperature step 
(<ta) do not get any better. We use circular patches for WMAP data analysis rather than square 
patches because of the computational advantage. Algebraic expressions for the observables given 
in equations (f6l)- (110p work very well for circular patches with negligible numerical differences for 
the case of square patches. A square patch with 28 x 28 pixels at the normal resolution of WMAP 
W-band data covers nearly the same area as the circular region with angular radius Or = 1.8°. We 
performed a further detecting power test with 28 x 28 pixels patches and found empirically that 
the relation between the critical value of A, A c , above which "good" results start to come out, and 
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a G is 

A c ~ 0.25o- G . (12) 

The different choices of input values for the parameters p, 6, or To showed no noticeable difference 
in the evolutions of errors. 



2. Application to WMAP 

An observed CMB anisotropy map is an aggregate of various independent modes of perturba- 
tion ranging from tiny sub-horizon scales to super-horizon scales well beyond the correlation length. 
As illustrated in Figure [IJ fluctuations with larger or smaller scales compared to the size of a test 
patch are neatly prescribed, but the intermediate scales whose wavelengths are comparable to the 
size of a test patch will appear to be continuous temperature tilts that also give plausible values 
of the step parameters. One drawback of this algorithm is that it does not distinguish between 
a discrete temperature step and a smooth temperature slope. However, this shortcoming can be 
easily fixed: If an apparent temperature step is detected at a spot on the map, we repeat the 
analysis with a half-sized patch at the same spot. If the structure is a continuous slope, then it 
would return half the value of A (, s than the previous result, while for the signal from a discrete 
step, the returned A t> s should stay the same within the error. We conducted the step signal search 
through the WMAP 3 year W-band data set and found 193,160 unqualified signals against back- 
ground fluctuation (the constraint in eq. [12] is not applied), 129,049 qualified steps+tilts (signals 
that meet the constraint eq. [12j : i.e., A. b s > Q.25aa), and 12,330 qualified discriminated steps. 
The A D ft s value of qualified steps is in the range 

(18.3 ± 3.0)fiK < A obs < (115.4 ± 6.0) fiK (13) 

as the curves show in Figure El Patches with radius Or = 1.8° in the WMAP 3 year W-band have 
a background fluctuation oq less than 176 /iK in the Kp2-mask-cleared region. This means, at 
the worst case, we can identify a step signal as low as 176 /uK/4 = 44 with moderate errors. 
Therefore, if there are cosmic string signals with A > 115.4 //K and if they are located in the 
available region (out of the Galactic plane or the Kp2 masked region), they would not be missed. 
Thus, we can set an upper limit on the cosmic string signal 

A = 4nGn<y 8 8 T\ cos <f>\ < 127 AfiK, 95% CL (14) 



where <j) is arbitrary. Thus, with ("fsPs) — 1 (jVilenkin &: Shellard 19941 ). the upper limit of the 
cosmic string parameter G\i can be estimated as 

Gfi < 3.7 x 10~ 6 , 95% CL. (15) 
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3. CONCLUSIONS 

We have presented and tested a new technique to directly discover cosmic strings via the 
patterns they would produce in a CMB anisotropy map. We found that the minimum magnitude 
of the step signal that is required to be unequivocally identified is 44 fiK for the WMAP 3 year 
W-band data set. This algorithm can be used to crop the reliable step signals from the CMB 
anisotropy data, and it will serve as the valuable ground work for future pattern searches with 
more refined data, such as further WMAP data releases or PLANCK data. 

Computer simu lations and data analysis with the WMAP data set were done using HEALPbfl 



(jGorski et al. 20051 ). This work was supported by LBNL and the Department of Physics at the 



University of California, Berkeley. 
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Fig. 2. — Evolution of errors of characteristic parameters as functions of A with input values of 
parameters (To = W.O^K, p = 0.5, 9 = 0, og = 116-OfiK). Left panels: 8x8 pixels in a square. 
Middle panels: 14 x 14 pixels. Right panels: 28 x 28 pixels. For each size of patch, the number of 
simulations performed is 10,000. The shaded regions indicate the ranges of A for which the step 
feature in a patch is not obvious compared to its background fluctuation (ac = 116.0 /xK). The 
result on A f, s is biased upward because the simulations with erroneous values of A a b s (A b s < 0) 
were not included in the statistics. 
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Fig. 3. — Distributions of A Q b s detected in WMAP 3 year W-band data set. The number of step 
signals found are 193,160 unqualified signals [dotted curve), 129,049 qualified steps+tilts {dashed 
curve), and 12,330 qualified steps {solid curve). 



